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Abstract 

The NA60 experiment has studied J/tjj production in Indium-Indium collisions at 
158 A-GeV. In this paper we present an updated set of results obtained with the com- 
plete set of available statistics and an improved alignment of the vertex tracker. The 
centrality dependence of the J/V> production, obtained with an analysis technique 
based only on the J/ip sample, indicates that a suppression beyond that induced 
by nuclear absorption is present in In-In collisions, setting in at ~80 participant 
nucleons. A first study of the systematic errors related with this measurement is 
discussed. We also present preliminary results on the J /tp azimuthal distributions. 
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1 Introduction 



The study of dilepton production in heavy-ion collisions at the CERN SPS 
provides some of the most interesting observables explored so far in the search 
for the quark-gluon plasma. In particular J /if) suppression is one of the most 
important signatures for the formation of a deconfined phase pQ. The NA38 
and NA50 experiments have extensively studied J/if> production in various 
colliding systems, including p-A, S-U and Pb-Pb [2J. Proton-nucleus data pro- 
vide an important reference, which describes the expected absorption for J /if) 
crossing cold nuclear matter. By comparing the centrality dependence of the 
J/if> yield observed in A- A collisions to this reference, one can look for anoma- 
lous J /if) suppression. In particular, NA50 has observed, in Pb-Pb collisions, 
that below a certain centrality threshold the J /if) production is well described 
invoking nuclear absorption as the only suppression mechanism; on the con- 
trary, above that threshold, an anomalous suppression sets in, and the J /if) 
yield becomes considerably lower than expected from the nuclear absorption. 
The observations done by NA38/NA50 surely give a clear indication for an 
anomalous J /if) suppression, however they should be complemented by other 
sets of accurate measurements obtained from different collision systems. The 
NA60 experiment at CERN SPS, with an innovative experimental apparatus 
improving the dilepton detection technique, studies J /if) production in In- In 
collisions, in order to investigate if the suppression is already present in sys- 
tems lighter than Pb-Pb. The comparison between different colliding systems 
should also allow to understand which is the physics variable at the origin of 
the J /if) suppression as a function of centrality. 

2 The experimental apparatus 

The NA60 experimental apparatus is based on the muon spectrometer inher- 
ited from the NA50 experiment, complemented by a new target region, based 
on a Beam Tracker and a Vertex Tracker placed in a 2.5 T dipole magnet. 
The target system is made of 7 In targets, 1.5 mm thick each, placed in vac- 
uum, corresponding to a total interaction probability of ~20%. A Zero Degree 
Calorimeter (ZDC) provides an estimate of the centrality of the collisions, by 
measuring the energy (E Z dc) released by the projectile nucleons which have 
not taken part to the interactions. A more detailed description of the appara- 
tus can be found in [3lH] . The reconstruction of the muon tracks in the muon 
spectrometer is affected by multiple scattering and energy loss fluctuations 
due to the crossing of the muons through the hadron absorber. To overcome 
this limitation, NA60 measures muons before their entrance into the hadron 
absorber, by means of the vertex tracker. Muon tracks measured in the muon 
spectrometer are matched with tracks reconstructed by the vertex tracker, re- 
sulting in an improvement of the dimuon mass resolution and in an accurate 
determination of the muons' origin. Details on the matching technique can be 
found in [5]. 
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3 J /if) results from the Indium-Indium run 

During the 5 weeks long In-In run in 2003, NA60 has collected ~230 million 
dimuon triggers, running with beam intensities around 5xl0 7 per 5-second 
spill. Two muon spectrometer settings have been used during the data taking, 
corresponding to two values of the magnet current: 4000A and 6500 A. Events 
collected with the lower current have a better acceptance for the low mass 
dimuons, while the higher current setting improves the mass resolution of the 
J /if). With respect to the previously published results [6], the analysis pre- 
sented in this paper includes both data sets, i.e. the whole statistics collected 
by the NA60 experiment. Furthermore, a new reconstruction of the data has 
been performed. It is based on a better alignment of the vertex tracker, result- 
ing in few fim accuracy, and on an improved quality of the tracking algorithm 
in the vertex spectrometer. Details on these issues can be found in [TJ. 

The event selection criteria have been tuned in order to obtain a clean sample 
of J /if) events produced in primary In-In collisions. Tracks are reconstructed in 
the vertex tracker and we require to have at least one vertex satisfying various 
quality cuts, including a minimum number of tracks (> 4) attached to it. The 
longitudinal coordinate of that vertex, which is reconstructed with an accuracy 
of ~200 /im, must lie inside one of the seven targets. The matching of the muon 
tracks makes it possible to determine the longitudinal coordinate of the dimuon 
production point with a precision of few hundreds /im. Asking the dimuon 
vertex to coincide with the most upstream interaction vertex in the target 
region ensures that the J /if) is produced in a primary In-In collision. A further 
set of selections is done in order to minimize the event pile-up in our apparatus. 
It is based on the good timing properties (a few ns) and on the excellent 
granularity of the Beam Tracker, that allows to reject events with several ions 
arriving close in time. This selection is particularly important for obtaining 
an unbiased Ezdc distribution. In fact, since the ZDC is placed on the beam 
axis and records a signal for each incoming ion, it is particularly sensitive to 
this kind of background. Finally, in order to reject events coming from the 
edges of the spectrometer acceptance we only retain dimuons corresponding 
to the phase space window < y cms < 1 and — 0.5 < cos 9c s < 0.5, where 9c s 
is the polar angle of the muons in the Collins-Soper reference frame. 

Two complementary analysis have been performed in order to study the cen- 
trality dependence of the J/if> production in Indium-Indium collisions. They 
differ in the way used for normalizing the J/if> distributions. The first tech- 
nique, the so-called "standard" analysis, is based on the normalization of 
the J /if) to the Drell-Yan (DY) yield. In this approach, already used by the 
NA38/NA50 experiments, Drell-Yan is used as a reference process. In fact, 
being independent from final state effects, its production cross section scales 
with the number of collisions and represents, therefore, an ideal normalization 
for a hard process, like 3 /if) production. In the "standard" analysis, the ratio 
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between the J/ip and the DY yields is extracted by fitting the opposite sign 
dimuon mass spectrum to a superposition of different contributions. Above 
2 GeV/c 2 , the dimuon invariant mass region contains the J/ip and the ip' 
resonances, an underlying continuum composed of Drell-Yan events and semi- 
muonic decays of D and D mesons, and a combinatorial background from 7r 
and K decays. The expected mass shapes of the signals are evaluated through 
a Monte Carlo simulation based on Pythia with GRV94LO parton distribu- 
tion functions. These simulations also provide the values of the J/i/j and DY 
acceptances (Aj/^—12A%, A DY (2.9-4 .5) = 13.0% for the data collected with the 
current of 6500A in the muon spectrometer and slightly larger for the 4000A). 
The combinatorial background has been estimated from the measured sample 
of like-sign pairs. Results obtained from the two sets of data, corresponding to 
the two configurations of the muon spectrometer, are found to be statistically 
compatible. Therefore in the following we show results based on their average. 
The study of the ratio between J ftp and DY cross-sections has the advantage 
of being free from systematic errors connected to experimental inefficiencies 
and to the integrated luminosity, since both J/ip and DY events were collected 
under the same conditions. However, as it was mentioned, the statistical error 
is large, due to the limited number of high mass DY. Therefore, in order to in- 
crease the statistics as much as possible, the matching of the muon candidates 
in the vertex spectrometer has not been required. This allows to increase the 
statistics by ~40%, for a total of ~320 DY events with mass above 4.2 GeV/c 2 
and -45000 J/ip. In this further cut based on the extrapolation of the 

muon spectrometer tracks ensures that the dimuon comes from the target re- 
gion. In spite of the statistics gain, the number of high mass Drell-Yan events 
is still limited; as a consequence the J/ip / DY ratio cannot be studied in more 
than three centrality bins. 




Fig. 1. Left: J /ip /DY standard analysis divided by the normal nuclear absorption 
curve vs. N par t. NA50 and NA38 results are also shown. Right: measured J/ip com- 
pared to the normal nuclear absorption curve vs. Ezdc- The error on the relative 
normalization between the absorption curve and the measured points is not shown. 
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The J /if} I DY values have to be compared to expectations from a pure nuclear 
absorption scenario (the so-called "nuclear absorption curve"). This reference 
distribution is obtained using the Glauber model, assuming that J/0 pro- 
duction is a hard process and adopting the value <r}^=4.18 ± 0.35 mb, as 
measured by NA50 [2J, for the nuclear absorption cross section of the cc state. 
In Fig. [1] (left) the measured J/0 /DY ratios have been divided out by such 
a curve and plotted as a function of the number of participants (N part ). This 
centrality estimator has been obtained from Ezdc using the Glauber model 
and taking into account the smearing induced by the detector resolution. We 
see that a suppression of the J/0 beyond nuclear absorption is present in In- In 
collisions. However, the small number of bins prevents an accurate study of 
the J /ifj centrality dependence as well as a meaningful comparison with the 
results obtained by NA38/NA50. 

The second analysis technique overcomes the problem of the low Drell-Yan 
statistics by directly studying the measured J/0 centrality distribution as a 
function of Ezdc- Following this approach, the J/0 sample is not normal- 
ized to a measured quantity. Rather, the observed J/0 centrality distribution 
is directly compared to the theoretical distribution expected in case nuclear 
absorption is the only suppression mechanism. Since the number of J /if} is 
not normalized to other data collected under the same conditions, one has to 
check that no centrality-dependent inefficiencies are introduced by the selec- 
tion criteria and, more generally, that systematic errors are reduced as much 
as possible. For this reason in this analysis we require the matching of the 
muon tracks in order to identify and discard events where the J/if> is not 
produced in primary In-In collisions. We have checked, by means of a Monte 
Carlo simulation, that the matching efficiency as a function of centrality for 
J /if} events increases by less than 2% from central to peripheral events, in- 
ducing a negligible bias on the centrality distribution of the measured J/ip. In 
addition, the matching improves significantly the J /if} mass resolution, from 
105 to 70 MeV, and reduces the combinatorial background from 3 to 1% in the 
corresponding mass region. We have also checked that the vertexing efficiency 
for events where a J/0 has been produced, decreases by less than 1% from 
central to peripheral collisions. After event selection, the J /if} distribution as 
a function of centrality has been obtained, in 1 TeV Ezdc bins, by means of a 
simple fitting procedure that allows to subtract the small amount of Drell-Yan 
and combinatorial background under the resonance peak. The J/0 centrality 
distribution is then compared to expectations from pure nuclear absorption. 
The relative normalization between data and the reference curve is not fixed a 
priori; therefore we assume the ratio between data and the nuclear absorption 
curve, integrated over centrality, to be the same as in the "standard" J/0 / DY 
analysis, i.e. 0.87 ± 0.05. 

The result is shown in Fig. [1] (right) as a function of Ezdc, while in Fig. [2J 
(left) we plot the ratio between the measured and the expected J/0 yield. 
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Fig. 2. Left: ratio between measured J /if) and the absorption curve vs. Ezdc (stars: 
standard analysis). The three most central bins of Fig. Q] (right) have been grouped 
together, to reduce statistical fluctuations. Right: ratio between measured J /if) and 
the absorption curve as a function of ~N part for S-U, Pb-Pb and In-ln collisions. 

It must be noted that events where a heavy nuclear fragment produced in a 
primary In-ln collision reinteracts in a downstream target have a biased Ezdc 
value. A Monte-Carlo correction (smaller than 2%) has been introduced to 
correct for the apparent shift towards central events induced by this effect. 
This analysis confirms the presence of an anomalous J /if) suppression in In-ln 
collisions, already seen in the "standard" approach. The departure from the 
expected normal nuclear absorption occurs at ~12 TeV and is followed by a 
rather flat behaviour for more central events. In Fig. [2] (right) the In-ln re- 
sults, rebinned in order to reduce statistical fluctuations, are plotted together 
with the previous results of NA38 /NA50 as a function of the number of par- 
ticipant nucleons. Although this comparison must be considered preliminary, 
due to the possible systematics connected with the use of different centrality 
estimators (transverse energy for NA38/NA50, Ezdc for NA60), a fair agree- 
ment between the various systems can be observed. As a side remark, it can 
be noticed that the In-ln suppression pattern shown in Fig. [2] (right) is com- 
patible with the occurence of a step-like suppression as a function of N part , 
once the smearing due to the E Z dc resolution is taken into account. Such a 
kind of study shows that the departure from the normal nuclear absorption 
behaviour occurs at N pari ~80. As can be seen in Fig. [2] this approach provides 
results with negligible statistical errors (< 2%). Concerning systematic errors, 
we have verified that their most important source is connected with the deter- 
mination of the nuclear absorption curve and with the determination of N part 
starting from Ezdc- The calculation of the absorption curve is based on the 
Glauber model. Therefore we have investigated the influence of the choice of 
the inputs on the result. Different Indium nuclear density distributions [8] with 
respect to the default parameterization [9] have been tested. This choice does 
not affect significantly the shape of the absorption curve, since the effect is 
smaller than 1% for Ezdc> 3 TeV. For very central events the effect is larger 
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and reaches ~15%. Concerning the centrality determination, the zero degree 
energy is well correlated with N part . However, the ZDC does not measure only 
spectator nucleons but also a small amount of energy released by the forward 
secondary particles, emitted in the calorimeter's acceptance {rj > 6.3). The un- 
certainty on this contribution, important only for very central collisions, has 
been conservatively assumed to be of the order of 10%. We have calculated 
that such an uncertainty gives, for events with Ezdc < 3 TeV, a systematic 
error of the order of 9% on the absorption curve. For more peripheral events 
the effect is negligible. Furthermore, the calculation of the absorption curve 
depends on cr^ and on both estimated in p-A collisions at 450 GeV by 
NA50. The experimental errors on such values, and in particular the factors 
needed to rescale a p fj^ from 450 to 158 GeV, induce a systematic error on the 
absorption curve of the order of 8%, almost independent of centrality. Finally, 
since we have normalized the data to the absorption curve using the results 
from the "standard" analysis, the 6% error related to the uncertainty in the 
normalization factor has also to be taken into account. Combining the various 
sources, we end up with a ~10% systematic error, independent of centrality. 
On top of that, the most central bin {Ezdc < 3 TeV) is affected by a fur- 
ther, sizeable systematic error relatively to the others. This means that the 
shape of the suppression pattern has been determined with good precision, 
while the absolute normalization of the points is more uncertain. It must be 
noted, anyway, that a significant fraction of the systematic errors originates 
from the uncertainty in the energy rescaling of the absorption curve. The p-A 
data collected by NA60 at 158 GeV should help to reduce significantly this 
uncertainty. Finally, one should be aware that the systematic errors related 
to the energy rescaling of the absorption curve also affect the NA38/NA50 
points. Therefore they do not influence the relative positions of the patterns 
observed in S-U, In-In and Pb-Pb. 

NA60 has also obtained preliminary azimuthal angle distributions of the J/ip 
produced in In-In collisions. These distributions have been obtained with 
respect to the orientation of the reaction plane, defined by the beam di- 
rection z and the impact parameter vector b. This plane is experimentally 
not directly accessible, but it can be determined from the emission angles 
of the charged particles measured in the vertex tracker, using the so-called 
event plane method. The J/ip anisotropic flow is then quantified evaluating 
the coefficients of a Fourier expansion of the azimuthal distribution. The 
azimuthal angle distribution for J/-0 is shown in Fig. [3] (left) for central 
events (0.5 < <r/a geo < 20%) and in Fig. [3] (right) for peripheral events 
(20 < cr/cr g eo < 83%). The elliptic flow coefficient v 2 is obtained after correct- 
ing for the event plane resolution. The distributions are based on ~12000 J/ip, 
which is about ~50% of the full statistics available for this kind of analysis. 
The more peripheral data seems to indicate a non-isotropic emission pattern 
but the limited statistics does not allow, for the moment, strong conclusions. 
A detailed review of this analysis is reported in [TP"] . 
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Fig. 3. J/ip azimuthal pattern for (left) central and (right) peripheral collisions. 
4 Conclusions 

The analysis of the data taken by NA60 in In-In collisions shows that the J /if) 
is anomalously suppressed. An analysis method based on the comparison of 
the J /if) centrality distribution with the expectation from a pure nuclear ab- 
sorption scenario shows that the anomalous suppression sets in at N part ~80, 
with a saturation for more central events. The statistical errors are negligible 
(of the order of 2%). We have estimated a 10% systematic error on this analy- 
sis, independent of centrality and therefore not affecting the shape of the 3 /if} 
suppression pattern. Most of the systematic errors come from the uncertainty 
in our knowledge of the J/if> nuclear absorption, which has been accurately 
measured at 450 GeV, but not at 158 GeV, the energy of the heavy-ion data 
taking. The situation is expected to improve as soon as results from the anal- 
ysis of the p-A data collected by NA60 at 158 GeV will be available. 
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